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Electromagnetic interactions are discussed in the context of the Klein- 
*/~) 1 Gordon fermion equation. The Mott scattering amplitude is derived in 

leading order perturbation theory and the result of the Dirac theory is 
reproduced except for an overall factor of sixteen. The discrepancy is not 
\ resolved as the study points into another direction. The vertex struc- 

O j tures involved in the scattering calculations indicate the relevance of a 

modified Klein-Gordon equation, which takes into account the number 
of polarization states of the considered quantum field. In this equation 
the d'Alembertian is acting on quaternion-like plane waves, which can be 
generalized to representations of arbitrary spin. The method provides the 
same relation between mass and spin that has been found previously by 
Majorana, Gelfand, and Yaglom in infinite spin theories. 
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1 Introduction 

Quaternionic calculation methods are under permanent investigation and fur- 
ther development to provide powerful mathematical tools for applications in 
various areas of physics. Silberstein pQ, Lanczos [2], Conway [3], Giirsey [HE], 
Finkelstein et. al. [6], Edmonds [7], Adler [8l|9l[T0], Horwitz and Biedenharn 
?H '. [H] , Baylis [T2J , and many other authors inspired with their research on this sub- 

ject other physicists. The number of publications concerning with quaternions 
or complexified quaternions is in fact endless. A bibliography with hundreds 
of entries has been compiled by Gsponer and Hurni in alphabetical |13j and 
analytical form [14J . Recent articles on quaternions in relativistic physics have 
been published by Demir et al. [13 [H], Panicaud [17], Bisht et al. [151 fT9]. 
Christianto et al. [2D], and Tamsh [2T] , 

In the following text a quaternionic Klein-Gordon equation will be applied to 
fcrmions. Some decades ago Feynman expressed his predilection for the Klein- 
Gordon fermion equation, because the representation of relativistic quantum 
mechanics with the Klein-Gordon equation by the method of path-integrals can 
be easily handled, but the Dirac representation is very hard to represent directly 



1 



[2"2"1 |2"31 [21] • The Klein-Gordon fermion equation, which is equivalent to the 
squared Dirac equation, has been derived by Kramers 25, 26^ and Lanczos |27j . 
Field theoretic aspects of this equation have been investigated by Brown [55] and 
Tonin [5U] . See also Pilkuhn [30. for a detailed historic overview. More examples 
of second order differential fermion equations can be found in Bicdcnharn et al. 
[3"T] . Volkovyskii [35], or recent articles published by Angeles and Napsuciale 
[53] , Nottale and Celerier [331 133 , and Delgado- Acosta and Napsuciale [35] . 

The Klein-Gordon method in [37] is based on paravectors, which include the 
Pauli algebra as basis elements. Congruent complexified quaternion representa- 
tions have been considered before also by Edmonds |7J, Baylis [38], Sobczyk [39] . 
and Sachs [40] [41] . Recently, the Pauli algebra has been applied to fermion spin 
by Baylis et al. [49 , to polarization optics by Tudor [55j , and to semiconductor 
physics by Dargys, see [5T] and references therein. 

Quaternions are used as building blocks in the representation theory of Clif- 
ford algebras. Starting point for the application of Clifford algebras in modern 
physics was the book of Hestenes [42] . For more details on further progress in 
this area see the recent articles of Hestenes [33] and Sobczyk [33] . More infor- 
mation on Clifford algebras in physics can be found in Doran and Lasenby [45) , 
Rodrigues and Capelas de Oliveira [35], Girard [37J, or Perwass [15] . 

In [37) quaternions have been complexified with the product of the complex 
and the hyperbolic unit to obtain the hyperbolic Pauli algebra. The hyperbolic 
numbers appear in fact in two representations. The complex number like rep- 
resentation has been considered by Shervatov [55] , Yaglom [53J [53] , Ryan [55J , 
Hucks [55] , Gal [57] , Yamaleev [55] [59] , Catoni et al. [55J [5T] , and many other 
authors. On the other hand there is the null plane representation in terms of the 
double field. An overview of the double field and its role in the representation 
theory of Clifford algebras has been provided by Porteous [62] . The hyperbolic 
numbers have been used by Khrennikov )63) and Nyman |64j to represent prob- 
abilistic data by quantum-like hyperbolic amplitudes. This research area has 
been detailed in [55J [55] . The work on the Erlangen program of Kisil considers 
complex, parabolic and hyperbolic numbers and their geometric properties, see 
for example [67j [68] [59] . Solutions of the Klein-Gordon equation with hyper- 
bolic numbers have been studied by Bracken et al. [70] and Kravchenko et al. 
[7T][7J]. 

The following article applies the Klein- Gordon fermion theory as presented 
in [37J [73] to calculate the Mott scattering amplitude. The result of the Dirac 
theory is reproduced except for a factor of sixteen. Though it might be possible 
to close the gap with a more detailed analysis, another method is proposed in 
Section [9J which is based on the structure of the interaction vertex. The Klein- 
Gordon equation is modified to include the number of possible polarization 
states of the considered quantum field. Quaternionic plane waves are introduced 
in Section [TU] as solutions of the equation, which could be denoted as Klein- 
Gordon equation for polarizable fields. 

The mathematical representation will be adjusted in Section [5J to align the 
formalism more closely with the Clifford algebra approach. Because the article 
should serve as a bridge between the method specified in [37] and the modified 
concepts, the notation is adjusted in the middle of the text and not from the 
beginning. 
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2 Hyperbolic Pauli Algebra 



The following considerations are based on the hyperbolic Pauli algebra, which 
is used to represent relativistic space-time vectors. A relativistic vector x is 
expanded into 

x = e^. (1) 

The basis elements are defined as a paravector algebra, which means one basis 
element of the algebra is the identity element. The other basis elements are 
constructed of the Pauli algebra multiplied by the hyperbolic unit j = y/+T 

= (e ,e k ) = (l,j<Tfc) . (2) 

The hyperbolic Pauli algebra e k is isomorphic to the Pauli algebra <j}~, but 
includes the hyperbolic unit as an extra factor. The multiplication rules of the 
hyperbolic Pauli algebra are 

e k ei = 5 kl l + ije k i m e m . (3) 

In the same way as the Pauli algebra can be understood as the complexification 
of the quaternions by the complex unit i, the hyperbolic Pauli algebra can be 
understood as a complexification of the quaternions by the factor ij. 

The hyperbolic Pauli algebra corresponds to the real Clifford algebra B^o- 
The advantage of the hyperbolic Pauli algebra, compared to the standard Pauli 
algebra, is that the Clifford anti-involutions conjugation and reversion and the 
grade involution can be defined in a convenient form, when using a 2 x 2 matrix 
representation of the algebra. Conjugation as matrix transposition and change 
of sign of the hyperbolic and complex units, reversion as matrix transposition 
and change of sign of the complex unit, and the grade involution as change of 
sign of the hyperbolic unit [73] . 

In the mentioned 2x2 matrix representation of R3,o the scalar product can 
be defined with the conjugation anti-involution as 

x ■ y = ^Tr(xy) . (4) 

Though the scalar product within Clifford algebras is normally defined in a 
representation free form, the trace definition will be used to get a deeper un- 
derstanding of an expression that appears in the scattering calculations of the 
following sections. 

3 Field equations and interactions 

The Lagrangian for free noninteracting spinor fields has been defined in |73] as 

C = 'tp(M 2 - to 2 )V> . (5) 

The mass operator is represented in terms of the relativistic paravector algebra 
introduced in Section [5] 

M 2 = pp . (6) 

The coefficients of the momentum operator are p^ = id^. They are invariant 
under conjugation 

= • ( 7 ) 
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The equation of motion, which is equivalent to the Klein-Gordon equation, can 
be derived from the above Lagrangian according to 

M 2 tp = m 2 ijj ■ (8) 

Interactions can be introduced as usual with the minimal substitution. The 
covariant derivative is defined as 

P=p + A, (9) 

which ensures the invariance under gauge transformations. A coupling constant 
has been omitted for simplicity. The Lagrangian in Eq. (|S|) remains formally 
invariant after the minimal substitution. However, the mass operator 

M 2 = PP (10) 

is now equivalent to the quadratic Dirac operator with interaction terms [37] . 

The field equation of the gauge boson can be represented with the mass 
operator as given in Eq. ^ , which is now acting on a vector field 

M 2 A = -J . (11) 

This representation includes all Maxwell equations [T^J EZ] • The current J is 
the source term of the interacting field. Contracted with a polarization vector 
e the current can be written as 

e ■ J — ip(p + A)eip — $s(p — A)ip . (12) 

Note that the momentum in the second term on the right side of the equation 
is acting to the left. Using the plane wave expansion of defined in [37], the 
matrix elements of the current can be deduced in leading order perturbation 
theory as 

(f\e- J\i) =u f (pie + epf)ui . (13) 

The applied Feynman rules remind of scalar electrodynamics, but they include 
spin matrices and spinors. 



4 The scattering amplitude 

In this section the scattering amplitude is defined and further investigated in 
the context of unpolarizcd particle beams. The scattering amplitude may be 
introduced in the following form 

e 

with the sum over the possible polarizations of the exchanged gauge boson. The 
indices indicate vertices 1 and 2. The propagator function of the gauge boson 
is denoted as Da(q)- The absolute square of the scattering amplitude may be 
written as 

\M\ 2 = \D A (q)\ 2 J2^2, (15) 
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where the contributions arising from the two currents have the structure 

AT= (e • J)(f • J) . (16) 

The indices have been dropped for clarity. If the spin polarization of the scat- 
tered particles is not considered, one has to average the spin contributions of the 
initial states and take the sum of the outgoing states. To calculate the squared 
scattering amplitude in first order perturbation theory without exchange term 
the current as given in Eq. (fT3|) is inserted 

(A/ - ) = ^ ^2 Ui{pf£ + epi)ufUf(pi^ + spf)ui . (17) 

Si Sf 

The sum of the final spin states results in the identity element due to the 
completeness relation of the spinors [37] 

W = 2 u ^Pf s + £ Pi)(PiS + 0f)ui ■ (18) 

Si 

For the initial states the Casimir trick can be used to simplify the expression 
according to 

(A?) = \Tr[(p } e + . (19) 

The calculations have been performed so far based on the polarization vectors 
of the gauge boson to provide a tensor free representation. This notation will 
not be applied further in the following text. 



5 Momentum dependence of the scattering am- 
plitude 

In order to make the following calculations easily accessible, they are performed 
in a conventional tensor representation. The squared absolute scattering ampli- 
tude is now written as 

\M\ 2 = \D A {q)\ 2 Nx^N$ v , (20) 

where the completeness relation of the polarization vectors has been taken into 
account. The involved tensors are defined as 

Mpu = J»Jv ■ (21) 

The tensor free matrix elements of Eq. (fT"3|) will be replaced by a current of the 
form 

(/I J n N) = Ufipie^ + e^p^Ui . (22) 

Before the calculations proceed, some useful relations are provided. In anal- 
ogy to the notation used by Doran and Lasenby |45j one may introduce the 
relationship 

Cfi^is — ' ^-v &[X A &u • (23) 
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0110 


1001 


0220 


2002 


0330 


3003 


-1 


1221 


2112 


2332 


3223 


1331 


3113 


ij 


0132 


3201 


0213 


1302 


0321 


2103 


ij 


1023 


2310 


2031 


3120 


3012 


1230 


-ij 


0123 


2301 


0231 


3102 


0312 


1203 


-ij 


1032 


3210 


2013 


1320 


3021 


2130 



Table 1: Non-zero entries of Tr (cr p/J aw) /2. The first column shows the non- 
zero value, the other columns the indices, which have this value 



The right side of this equation can be further transformed into a conventional 
tensor representation, which is useful for explicit calculations. The scalar prod- 
uct of the basis vectors is given as 



The wedge product is represented according to 

The spin tensor, including the Pauli matrices Cj, is calculated as [57] 
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(24) 
(25) 

(26) 



With the help of this notation, the tensors involved in the scattering ampli- 
tude are calculated for unpolarized particle beams in analogy to Eq. (|19[) 



(27) 



The spin sums have been evaluated using the Casimir trick and the completeness 
relation of the spinors. Initial and final momenta have been replaced by the 
following definitions 



d = Pf + Pi , 



Pf-Pi 



(28) 



Because two terms in Eq. (f2"7T) vanish due to the trace operator, the tensor can 
be represented in the form 

(A/V = d^d u - . (29) 
The second term on the right side of the equation is an abbreviation for 



(30) 



The trace has non-zero entries as displayed in Table [TJ The spin contributions 
thus provide the following structure 
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(31) 



G 



The non-zero entries of the tensor r\p. vpa correspond to the first two rows in 
Table[T] The second term is proportional to the anti-symmetric four-dimensional 
Levi-Civita symbol. Thus the tensor defined in Eq. (I3U1) can be calculated as 

= r i P ^aq p q a ■ (32) 

The antisymmetric contributions of the tensor vanish due to the symmetry of 
the two momentum vectors. 

Now Eq. (|29[) can be inserted into Eq. (|2T)]) and the squared absolute scat- 
tering amplitude of the one-boson exchange Feynman diagram is calculated as 

\{M)\ 2 = \D A (q)\ 2 [(di • d 2 ) 2 - {d x d x + d 2 d 2 ) -b + b-b}. (33) 

The dot product has been generalized to tensors of rank two. The notation for 
the mixed terms is resolved according to 

drd t ■ b = r] p ^qPd»d»q° . (34) 

As a final remark it should be mentioned that the tensor b can be written as 
the scalar product 

= ( e M A i) ■ ( e » A l) (35) 
where Eqs. (U]) and (|23|) can be used to resolve this notation. 



6 Mott scattering 

Based on the results of the previous section it is now possible to derive the 
squared absolute amplitude for an electron scattered of a much heavier particle, 
for example a proton. The proton with mass M is at rest in the lab frame. The 
energy of the electron E is related to the electron mass m. 









( M \ 




( E \ 

















Pit = 





P2i = 





Plf = 


p sin 9 




V p J 




V o j 




i pcos9 J 



(36) 



Energy and momentum of the outgoing proton are determined by energy-momentum 
conservation. 

With this parameterization one can calculate the vectors which contribute 
to the scattering amplitude of Eq. (|33p . Only terms proportional to the proton 
mass are considered in the following, because the proton is very heavy compared 
to the electron. Thus one finds 

(di • d 2 ) 2 w {AEMf = 16M 2 (m 2 + p 2 ) . (37) 

The second and the last term of Eq. (13U|) can be neglected, because these terms 
do not include any contributions proportional to the proton mass 

dxdi ■ b w b ■ b « . (38) 

In the third term only the energy coordinate of d 2 contains the proton mass. 
With Eq. (f34|) and the 77-tensor as given in the first two rows of Table [T] one 
finds 

d 2 d 2 ■ b « (2M ) V = 16MV sin 2 9/2 . (39) 
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Inserting now all contributions into Eq. (|3"3")l the following result is obtained 

\(M)\ 2 = \Q\D A {q)\ 2 M 2 {m 2 + p 2 cos 2 9) . (40) 

The scattering amplitude needs to be multiplied by a phase space factor to 
derive the Mott formula. However, the phase space calculation rules of the 
applied method have not been investigated so far. 



7 Dirac Theory 

The calculation of the Mott formula within the Dirac theory can be found in 
many textbooks, see for example Griffiths [7J]. The corresponding result for the 
scattering amplitude is 

\(M)\ 2 = \D A (q)\ 2 M 2 (m 2 +p 2 cos 2 9) . (41) 

The result differs by a factor of sixteen from the result of the Klein- Gordon 
fermion theory in Eq. (|4U|) . The discrepancy could be resolved by a more in 
depth analysis of the Klein-Gordon approach. Because the mass operator in 
Eq. (fTU)) is equivalent to the squared Dirac operator both methods should pro- 
vide the same cross sections. However, even a solution of the factor problem 
is not able to invalidate the argument that the Klein-Gordon fermion theory is 
only an alternative representation of the Dirac theory and does not provide any 
new or deeper insights into physics. 

The rest of this paper thus considers a modified method, which is motivated 
by the mathematical structures that appeared in the scattering calculations of 
the previous sections. These calculations further indicate that the method speci- 
fied in [37] should be aligned closer with other Clifford algebraic and quaternionic 
approaches, see for example Adler [TU] or Baylis |12) . 



8 Alignment with Clifford algebraic methods 



In the theory of Clifford algebras the pseudoscalar is obtained by multiplication 
of all available basis elements. The pseudoscalar may be introduced within the 
hyperbolic Pauli algebra representation of R 3i o as 



i = eie 2 e 3 . 

Using the basis elements as defined in Section [5] one finds 

% = -ij . 

Equation ([3]) can be rewritten in terms of the pseudoscalar 

e k ei = -Sul + ie k i m e m . 



(42) 
(43) 
(44) 



This equation indicates how the hyperbolic Pauli algebra is applied in its rel- 
ativistic context. The next modification is the redefinition of the tensor in 
Eq. (J2SD 



(45) 
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ie 3 


ie% 







-en 


ie 2 


S3 





ei 


ie 3 




-e\ 
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The elements of the Pauli algebra have been replaced in the spin tensor by 
ei = j<Ji. The definition of the wedge product for the basis vectors needs to be 
adjusted accordingly. 

e^he v = icr^ . (46) 
The generators of the spin angular momentum are introduced as in |37j 
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^ - T , (47) 



but now based on the tensor as defined in Eq. (|45|) . Within this represen- 
tation the Pauli-Lubanski vector may be defined as 



t w <j V s pcr p v = la^up" = e M A p . (48) 



Note, that the Pauli-Lubanski vector is multiplied essentially by a factor of 
two compared to definitions in common textbooks like Itzykson and Zuber [75j . 
The last expression in Eq. (|48|) represents the coordinates of the Pauli-Lubanski 
vector as biparavectors, which can be considered as planes formed by the mo- 
mentum and the basis vectors of the coordinate system. See Baylis |12j for more 
details on the geometric interpretation of biparavectors. 



9 Klein-Gordon equation for polarizable fields 

The following considerations are based on the observation that the Pauli-Lubanski 
vector is hidden in the vertices of the previously investigated scattering prob- 
lem. One may have a look at the matrix elements of Eq. (fTB")) . For example 
the second term on the right side of the equation gives rise to the following 
representation 

ep = e ■ (p + w) , (49) 

which can be derived with Eq. ((23)) and the above definition of the Pauli- 
Lubanski vector w. As the combination of momentum and Pauli-Lubanski 
vector appears in the current, one could ask whether this structure plays a 
more important role in the momentum representation of the wave equation. 
Thus instead of using the mass operator in its original form, one could think of 
an operator defined as 

O = + w) ■ (p + w) . (50) 

One has to keep in mind that in this notation a conjugation is hidden, which 
evaluates according to 

w M = -Wfj, ■ (51) 

Together with the orthogonality of Pauli-Lubanski vector and momentum the 
following straightforward calculation is leading to 

O = p 2 - w 2 = m 2 + Am 2 s(s + 1) = (win) 2 , (52) 

with n the number of possible polarization states of the considered quantum 
field 

n = 2s + 1 . (53) 

Equation (|5"2"|) is valid for the n = 2 spin representation, but it can be generalized 
also to an arbitrary number of polarization states, including n = 1 for spinless 
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particles. The wave equation in a momentum space representation could have 
the compact form 

Oip = (nmfip . (54) 

This representation assigns the matrix structures in the wave equation to the 
Pauli-Lubanski vector. 

To give a representation of the wave equation in position space, the d'Alcmbcrtian 
is introduced as 

□ = d ■ d . (55) 
The Klein-Gordon equation for polarizable fields is then defined accordingly 

[□+H 2 ]i/) = 0. (56) 

For spinless particles the equation reduces to the ordinary Klein-Gordon equa- 
tion. The plane wave solution ip will be discussed in the following section. 

One may compare the modified wave operator with the mass operator intro- 
duced in Eq. ((6]). In the new notation one finds 

M 2 =p-(p + w) . (57) 

Without interactions the spin structure collapses to M 2 — p 2 due to the or- 
thogonality of momentum and Pauli-Lubanski vector. After interactions have 
been introduced one obtains the squared Dirac operator with interaction terms 
[37) . Thus one may understand Eq. (|57[) in a wider sense as a synonym for the 
Dirac method. By comparison with Eq. (IBTfl) the difference to the modified wave 
operator becomes obvious. 

The difference between the approaches can be explained alternatively if one 
traces the generalized mass term back to the initial ingredients and writes 

- n(E 4 ) = C 2 {E 4 ) =p 2 -w 2 = C 2 (T 4 ) + C 2 {E 4 /T 4 ) . (58) 

In this equation the quadratic Casimir operators with respect to the Poincare 
group £4, the translational group T 4 , and the coset space E4/T4 are displayed. 
Equation (l58l) has been introduced in analogy to the Laplacian used by Meng 
[76] in investigations of the Quantum Hall Effect in higher dimensions. To 
make the analogy with Meng more obvious one may consider the d'Alembertian 
of the ordinary Klein-Gordon equation as a Laplacian acting on the group of 
translations T4 in Minkowski space-time 

A (T 4 ) = -D(T 4 ) - C 2 (E 4 ) - C 2 (E i /T 4 ) . (59) 

This equation has the structure as discussed by Meng [75]. If one builds a 
mathematical method on top of a Laplacian in T4, but the real physical world 
in vacuum is £4, there is the risk that additional mathematical structures have 
to be introduced at a later stage to align the T4 method with experimental 
observations. This gives an argument to redevelop the methodology with the 
E4 Laplacian of Eq. (|58p from scratch. 

10 Quaternion waves 

Adler uses the notion of quaternion waves |77j which should be adopted in 
the following text for the solutions of Eq. (|56p . The term quaternion must be 
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understood here in a generalized sense. In fact it would be better to speak of 
quaternion-like plane waves. They are introduced as 

ip — exp (—ixp) . (60) 

This structure can be transformed into various mathematically equivalent rep- 
resentations. A representation referring to the Clifford methodology can be 
derived with the help of Eq. ([2"3"]l 

V> = exp {—i{x ■ p + 1)) , (61) 

with the Clifford angular momentum [45] 

I — x A p — — (xp — px) . (62) 

For the representation of spinless particles the relativistic angular momentum 
term can be discarded. 

Another view on the plane wave is given by a tensor representation of the 
form 

tf> = exp (-KBjrfj" + s^F") , (63) 
with the relativistic orbital angular momentum 

= jjy _ x » p n ( 64 ) 

and the spin angular momentum as defined in Eq. (T4"T)) . Compared to the 
standard plane wave, which is related to translations only, the quaternion wave 
includes translations, rotations, and boosts. The representation which is most 
useful to reproduce the momentum space operator of Eq. (|50|) is 

ip = exp ( — ix ■ (p + w)) , (65) 

where the d'Alembertian of Eq. (|55p is applied to the quaternion wave. 

The Pauli-Lubanski vector does not provide any new information once the 
momentum and the number of possible polarization states of the considered 
quantum field are given. The proposed method thus describes states of plane 
quaternion waves by the mass, the number of polarization states, and the mo- 
mentum four-vector 

\ip) = \m,n,p) . (66) 

The number of polarization states may be interpreted as the intrinsic dimension- 
ality of the quantum field. The quaternion waves refer to the Poincare group as 
the gauge group of the method. Spinors can be introduced in this context as el- 
ements of the vector bundle, which is associated to the Poincare group principal 
bundle. 

The exponent in Eq. ((60)) can be transformed to a different overall sign, to the 
opposite order of momentum and position vectors, or to the opposite location of 
the bar symbol. The second transformation changes the sign of the relativistic 
wedge product, the third transformation the sign of a product of the form 

x I p = — (xp — xp) . (67) 

The eight possible discrete states of the quaternion wave must be encoded in 
Eq. (l66l) when the mentioned three types of discrete transformations are con- 
sidered. The transformations form the cyclic group Z2 ® 1i <8> Z2, which is 
isomorphic to the CPT group. 
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11 Higher spin 



The combination of Eq. (1531) with the constant term in Eq. (|56[) provides the 
same relationship between spin and mass as initially found by Majorana |78j 
and later also by Gelfand and Yaglom [79] 



See here also the work of Fradkin [50] , Casalbuoni [5T] , and Bekaert et al. [52] ■ 
The constant k has been introduced as 



Equation ((68)) points towards an equivalence of the Klein-Gordon equation for 
polarizable fields with the mentioned methods, which also refer to an arbitrary- 
component relativistic wave equation in comparison to the finite-component 
Dirac equation. 

Relativistic wave equations for particles of arbitrary spin have been intro- 
duced also by Dirac [83], Fierz, and Pauli [84] [85]. Higher spin quaternion 
waves have been investigated by Edmonds 86 . Varlamov considers higher spin 
in combination with the CPT group [57] . 

12 Summary 

The Klein-Gordon fermion theory in the mass operator representation has been 
applied to electron-proton scattering. It has been shown that the scattering 
amplitude calculated with this method is in agreement with the Dirac theory 
except for an overall factor of sixteen. In principle, the results provided by the 
mass operator should be identical due to the equivalence of mass and squared 
Dirac operator. Further investigations are necessary to explain this discrepancy. 

The structure of the vertex, as provided by the analysis of the scattering prob- 
lem, indicates the relevance of a modified Klein-Gordon equation for polarizable 
fields. The mass is multiplied in this equation by the number of possible po- 
larization states. The modified Klein-Gordon equation cannot be aligned with 
a squared Dirac operator, but equivalence to concepts derived by Majorana, 
Gelfand, and Yaglom is indicated. 
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